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A quantum system with finite degrees of freedom can simulate a composite of a system 
and its environment if the state of the hypothetical environment is randomized by external 
manipulation. We propose to examine various techniques in quantum information processing 
in virtual noisy environment with two concrete examples. One simulates phase decoherence 
of a single qubit in a transmission line, and the other does in a quantum memory. In both 
cases, the bang-bang control, a typical example of useful techniques employed in quantum 
information processing, is observed to be effective to suppress decoherence. 
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1. Introduction 

Quantum computing currently attracts a lot of attention since it is expected to solve some 
of computationally hard problems for a conventional digital computer. 1 Numerous realizations 
of a quantum computer have been proposed to date. Among others, a liquid-state nuclear 
magnetic resonance (NMR) quantum computer is regarded as most successful. Demonstration 
of Shor's factorization algorithm with NMR is one of the most remarkable achievements. 2 
Liquid state NMR will be denoted simply as NMR throughout this paper. 

Although the current NMR quantum computer is suspected not to be a true quantum 
computer because of its poor spin polarization at room temperature, 3 it still works as a 
test bench for a more realistic quantum computer. For example, we have demonstrated our 
time-optimal implementation of two-qubit quantum algorithms by using an NMR quantum 
computer. 4,5 

A molecule employed in NMR experiments can be arranged to work not only as a quantum 
register but also as a composite system of a quantum system and its environment. It is 
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possible to introduce decoherence phenomena in the quantum subsystem by manipulating 
the environment subsystem. Moreover the combined system can be employed as a test bench 
to develop techniques to protect a quantum system from decoherence. Decoherence is one of 
the primary obstacles in constructing a working quantum computer and must be suppressed 
somehow. Decoherence effect in a quantum register has also been studied elsewhere. 6 

The main purpose of this paper is two-fold. Firstly, we show that decoherence can be 
generated by manipulating the artificial environment. Secondly, we verify by NMR experi- 
ments that decoherence control methods, such as a bang-bang control, 7,8 actually suppress 
decoherence. It should be noted that demonstration of the effectiveness of decoherence control 
methods is rather difficult in other systems due to their extremely short coherence times. 

Section 2 is a brief review of the theory of a quantum channel, which is a useful formalism to 
describe decoherence in a general context. In Section 3 we describe decoherence of a one-qubit 
system in terms of a quantum channel. There we discuss the method to suppress decoherence 
by a bang-bang control. In Section 4 we show that a two-qubit system may be regarded as 
a composite of a system (qubit 1) and an environment (qubit 2). We introduce an artificial 
environment by manipulating qubit 2, which causes decoherence in qubit 1. We show two 
illustrating examples and calculate decoherence rates in these environments. We also show 
that application of a bang-bang control to qubit 1 suppresses decoherence in both cases. In 
Section 5 we report the results of our experiments, which support our theory. Using a two-spin 
molecule we demonstrate the generation of decoherence and its suppression by the bang-bang 
control. Section 6 is devoted to summary and conclusions. 

2. Decoherence 

2.1 Quantum Channel 

Decoherence is an irreversible change of a state of a quantum system which has quantum 
correlation with its environment. The change of the state of the system becomes irreversible 
due to our lack of knowledge about the state of the environment. 

Decoherence is formulated in terms of a channel or a quantum operation 1, 9-11 as follows. 
Let M's and J4? e be the Hilbert spaces of the system and the environment, respectively. The 
initial state of the system is represented by the density matrix p s while that of the environment 
by p e . The state of the whole system changes following the time-evolution law, 

p s ® p e -> Up s ® PcU 1 . (1) 

Here U is a unitary operator acting on the Hilbert space of the composite system <g> ,3%,. 
We consider the case in which the initial state is a separable state p = p s ® p c . The states of 
the system and the environment are correlated via the transformation (1). Needless to say, 
the unitary transformation (1) is a reversible change. If we are interested only in the state 
of the system, the measurement outcomes are completely described by the reduced density 
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matrix 

p' s = £(p s )=Tr e (Up s ®p e U^, (2) 

where the symbol Tr e denotes the partial trace over J%>. The partial trace operation is non- 
invertible and the associated loss of information leads to decoherence. Even if the initial state 
p s is a pure state, the transformed state $(p s ) becomes a mixed state in general. 

The mapping p s — > <§(p s ) is called a channel or a quantum operation. 1 A channel is the 
most general mathematical device to describe changes of a quantum state, including unitary 
time-evolution, measurement process, decoherence and so on. It is known that for a channel 
there is a set of operators {E^} acting on Jf s such that 

g(p s ) = Y,E kPs El (3) 

k 

J2 E k E k = Is- (4) 

k 

Here I s is the identity operator on J%. Equation (3) is called an operator-sum representation 
of the channel S. Equation (4) implies Tr s &(p s ) = Tr s p s , Tr s being the trace over J^j, and 
hence it is called the trace-preserving condition. 
2.2 Mixing process as a quantum channel 

There is another approach to defining channels without resort to partial trace over the 
Hilbert space of environment. Assume that we have a set of unitary operators {Uk}, which 
act on and that we have a set of real numbers {pk} such that < pk < 1 and YlkPk = 1- 
We then define a transformation of the system density matrix p s by 

p s -> JZ(ps) = ^PkUkPsUl- (5) 
k 

They satisfy the condition (4) if we put = yjpu Uk ■ This argument tells us that if we apply 
a set of time-evolution unitary operators {Uk} on the system with a probability distribution 
{pk}, we will observe a decoherence- like phenomenon after taking an average of the measured 
data over k. We call the transformation M a mixing process. 

Although a mixing process is defined superficially without referring to an environment, it 
is mathematically a special case of a channel that is defined through interaction between a 
system and an environment. For given sets of unitary operators {Uk} and probabilities {pk}, 
we can construct a Hilbert space Jt? e = {Y2k c k\k}} by demanding formally that {\k)} is a 
complete orthonormal set. Moreover, we define an environment density matrix 



p 



c 



J2 Pk \k){k\ (6) 



and define a unitary operator U = J2 k Uk®\k)(k\ that acts on J^ s (g) J^ e as 

u(\A)®\k)) = (Uk\A))®\k). (7) 
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By substituting them into the defining equation of a channel (2), we obtain the mixing process 
(5). In this paper we use the mixing process as a procedure to build a channel. 

3. Decoherence in one-qubit system 

3.1 Phase flip channel 

Here we introduce the phase flip channel, which is a typical example of a channel. We take 
a one-qubit system and a one-qubit environment for simplicity. Assume that the initial state 
of the environment is 



IV-e) = VP|0) + y/T=p\\) (8) 

with a real number p (0 < p < 1). Here the vectors {|0), |1)} are eigenstates of o z such that 
a z \0) = |0) and cr z \l) = — 11), where a x ,y,z are conventional Pauli matrices. We take a unitary 
operator 

(\ \ 
10 
10 
\0 -1J 

which acts on C 2 <8> C 2 , where I is the two-dimensional identity matrix. By substituting them 
into Eq. (2) we obtain a channel 

S(p s ) = Tr e (u Ps ® |Ve)<^e| ^) = E oPs El + Exps^ 

with 



[7 = 7® |0)(0| +a z ®\l){l\ = 



(9) 



(10) 



Eo = (0\U\il> e ) = y/pI, 
E 1 = (l\U\iJ; e ) = ^l~pa z . 
Any initial state of a one-qubit system is parametrized as 



Poo Poi 



-(/ + a^ffj; + a^cr^ + a z a z ) 



(11) 
(12) 

(13) 



with real numbers a x ,a y ,a z such that a 2 + a 2 + a 2 < 1. The vector (a x ,a y ,a z ) is called the 
Bloch vector. The angle (j> defined by 

2pio = a x + ia y = e l ^\J~al 



i 2 + a 2 



(14) 



denotes the azimuthal angle of the Bloch vector, and is called the phase of the spin. The 
complex quantity a x + ia y is called an amplitude in the context of NMR. The channel (10) 
transforms p s to 



&(p B ) = ^p(I + a x a x + a y a y + a z a z ) + ^(1 - p)(I - a x a x - a y a y + a z a z 

( Poo (2p - l)poA 

l(2p-l)pi p u I 



(15) 
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The first expression in the right hand side shows that the phase is left unchanged with prob- 
ability p while it is flipped as e 1 ^ — ► e*^ +7r ) = — e 1 ^ with probability 1 — p. Hence it is natural 
to call this a phase flip channel. In particular, whenp = ^, the transverse components (a x , a y ) 
of the Bloch vector, or the off-diagonal elements poi = Pic vanish after the channel is applied 
and the information about the phase is completely lost. However, the diagonal elements poo 
and pn, which represent populations of spins in the states |0) and |1), respectively, remain 
unchanged. Due to these properties, decoherence generated via the phase flip channel is called 
phase decoherence. 

It should be noted that different sets of an initial state of the environment and a unitary 
operator of the whole system may yield the same channel. Instead of Eq. (8) and (9), we may 
take a mixed environment state 

Pe =p\ + ){ + \ + {l-p)\-){-\= l -( ^ 2P ~\ (16) 

where |+) and |— ) are the normalized eigenvectors of a x with eigenvalues 1 and —1 respectively, 
and the controlled NOT gate 

(\ o\ 
10 
1 

\o o i oy 

Substituting them into Eq. (2) we obtain again the phase flip channel Eq. (15). 
3.2 Random fluctuating field 

It is possible to reproduce the phase flip channel without resort to the partial trace. This 
is done by introducing the mixing process defined previously. 

Let us consider a single spin Hamiltonian 

ff rf = -w(t)y, (18) 

where co(t) is a randomly fluctuating field. We have taken the natural unit H = 1. The time- 
evolution operator associated with the Hamiltonian (18) is the phase shift gate 

S(6) = e i9a ^ 2 (19) 

with 

0= I uj{r)dT. (20) 
Jo 

The phase 8 integrates the effect of uj{t) in the interval [0, t]. In the context of NMR, the 
phase shift gate is implemented with a longitudinal magnetic field or a scalar coupling with 
another spins as we discuss later. The phase shift gate acts on the density matrix as 



sw*si«n- r em . (2D 
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Given a probability distribution p(6) which characterizes the random fluctuating field, the 
mixing process is evaluated as 12 



Therefore the absolute value of the off-diagonal elements of ^(p s ) become smaller than those 
of p s . When the average (e~ ie ) is a real number, the map ^#(p s ) reproduces the phase flip 
channel (15). This applies when p(9) = p(—0) for example. 

We can further simplify the model without losing the essence of the random fluctuat- 
ing field model. Suppose that uj{t) takes only two values, ujq ± 5u>, with the corresponding 
probabilities p(±5u)). We simulate phase decoherence phenomena according to this simplified 
random fluctuating field model. 

Relaxation phenomena for nuclear spin qubits, including phase decoherence, have long 
been studied in somewhat different context in non-equilibrium statistical physics and em- 
ployed as probes in condensed matter physics. 13 Studies on a relaxation process induced by 
environmental fluctuations have played an essential role when we want to extract information 
on the environments. Conversely, we can evaluate the time evolution of a system provided 
that the structure of the environment as well as the interaction between the system and 
the environment are given. Spin relaxation has often been analyzed by the phenomenological 
Bloch equations. The equations are characterized by two parameters called the longitudinal 
relaxation time T\ and the transverse relaxation time T<i- The former characterizes energy 
relaxation whereas the latter describes phase decoherence. Relaxation phenomena can also be 
described by the method of master equation. This is called a Redfield equation in the field 
of spin relaxation. Both the Bloch and the Redfield equations are valid only for the so-called 
narrowing limit where the characteristic time r e of the environment is very short resulting in 
an exponential decay in the relevant spin variables. In contrast, the random fluctuating field 
model, though phenomenological, is applicable to any time scale ranging from the narrowing 
limit (t c — ► 0) to the slow modulation with large r e . Therefore, it is reasonable to take the 
random fluctuating field model as a basis for simulating phase decoherence phenomena. 
3.3 Suppressing decoherence by the bang-bang control 

Several groups 7 ' 8 have proposed and analyzed a useful technique to suppress decoher- 
ence, which is called a quantum bang-bang control. We briefly explain the principle of the 
bang-bang control. When a qubit system evolves in time, the interaction with its environment 
usually causes decoherence in the qubit state. If, however, time-evolution of the qubit could 
be reversed by some methods, the qubit returns to its initial state and decoherence would be 




(22) 



For any probability distribution p(6), 
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eliminated. Concerning the phase decoherence, a time-reversal operation can be simply im- 
plemented with a pair of -zr-pulses assuming that the state of environment remains unchanged 
between two pulses. The 7r-pulses around the x- and — x-axes transform the qubit state with 
unitary operators 

respectively. The phase shift operator S(9) = e ie<7z ^ 2 has the property 

V*S(P)V = S(-9) = S(9)-\ (24) 

Hence, by inserting a pair of 7r-pulses in a product of phase shift operators S(9) we get 

V^S(9)VS(9) = I. (25) 

Therefore, the state p s of the qubit comes back to the initial one as 

(V^S(6)VS(e))p s (V^S(e)VS(6)y = p s . (26) 

The phase shift is also canceled if a pair of 7r-pulses is inserted as 

S(O 2 )V*S(02 + 9 1 )VS{9 1 ) = I (27) 

showing that the locations of 7r-pulse insertions may be chosen rather arbitrarily. 

The time-evolution operator S{ojt) affects the qubit state if the interaction with environ- 
ment causes a phase shift proportional to time. Here a; is a parameter which characterizes the 
environment state and strength of interaction between the system and the environment. Let 
us introduce a time interval i b and put 2ntb = t with a positive integer n. Then we have 

S(iot h ) • S(wtb) • • • S(wtb) • S(wtb) = S(ut). (28) 

V v ' 

2n 

By inserting 7r-pulses we recover the initial state since 

V^S(ut b )VS{ut h ) • • • V^S(ut h )VS{ojt h ) = I. (29) 

This argument indicates that phase decoherence is suppressed by applying a regular sequence 
of 7r-pulses on the qubit. 

In a general circumstance, the environment state is not stationary and the phase shift is 
not proportional to time. Then S(ut) is replaced with 

U(t;to) = S u{r)d^j . (30) 

Even in such a general case, if u(t) dose not vary rapidly and remains almost constant during 
the short time interval t^, it is legitimate to use approximation 

U(t + 2t h ;t + t h ) «E/(t + t b ;*o) (31) 

and hence 

V ] U{t G + 2t b ; t + t h )VU{t + t h ; t ) w /. (32) 
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Therefore the phase shift will be mostly canceled by inserting 7r-pulses with a short interval t\, 
in the time-evolution operator (30) and hence the associated decoherence will be suppressed. 

4. Two-qubit system as a composite of system and environment 

4-1 Artificial environment 

Any system in an environment has a Hamiltonian of the form 



where H s and H e govern intrinsic behaviors of the system and the environment, respectively, 
while H se represents interaction between them. The relation of the system and the environment 
is schematically depicted in Fig. 1. Zurek, 14 for example, discussed a simple model where a 
one-qubit system is coupled to an n-qubit environment through interaction of the form <j z ®<j z . 

Suppose the interaction H se is so weak that its effect on the system qubit is negligible 
compared with that of H s for a certain time scale r. Assume further that the system consists 
of two subsystems, which are referred to as subsystems 1 and 2. Then the system Hamiltonian 
H s is decomposed as 



Here H\ and H2 govern intrinsic behaviors of the subsystems, while H12 describes interaction 
between them. Under this decomposition, we may regard the subsystem 1 as a new system and 
the subsystem 2 as an artificial environment. The subsystem 1 will exhibit a decoherence-like 
behavior if the subsystem 2 simulates an environment that has many degrees of freedom. 

Particularly, nuclear spins used in NMR have long relaxation times of the order of r ~ 10 s. 
Thus nuclear spins are almost isolated from the environment in the time scale r. In such a 
circumstance it is legitimate to regard some of the spins as an artificial environment for the 
other spins. 

Zhang et a/. 15 experimentally studied the spin dynamics of 13 C-labeled trichloroethane, 
which has three spins in a molecule. They regarded three spins as a composite of a two-qubit 
system and a one-qubit environment. They claimed that they observed decoherence in the 



H t — H s + H c + Hi 



(33) 



H s — H 1 + H2 + H\2- 



(34) 




Fig. 1. System and its environment. The system consists of two subsystems, 1 and 2. 
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two-qubit system. However, the artificial environment must have a large number of degrees 
of freedom to introduce irreversible decoherence-like behavior in the system while the one- 
qubit environment employed in their experiments had not enough degrees of freedom. Hence 
their system exhibited a periodic behavior and failed to introduce an irreversible change in 
the system. Ryan et a/. 16 overcame this difficulty by employing a seven-qubit molecule for 
simulating complex dynamics of an environment. 

Teklemariam et al. 17 proposed to apply a stochastic classical field 18 on the artificial envi- 
ronment to generate artificial decoherence. Although the artificial environment has only a few 
degrees of freedom, it can simulate an open environment if it is randomized by an external 
stochastic field. With this idea, they observed decoherence-like behavior by using NMR. Al- 
though we closely follow Teklemariam et al, we simplify the principle for generating artificial 
phase decoherence phenomena to the limit. Thanks to this simplification, we are able to ana- 
lyze the phenomena without numerical calculations and compare theoretical and experimental 
results directly. 

Several remarks are in order. An experimental work that may be called demonstration of 
"quantum bang-bang control" was only reported by Morton et a/. 19 to the best of our knowl- 
edge, where the time development of Rabi oscillation was undone by successive 7r-pulses. Their 
experiments are closely related with the "quantum Zeno effect" experiment by Itano et al., 20 
which was in turn based on the proposal by Cook. 21 Cook pointed out that demonstration of 
"quantum Zeno effect" is difficult and thus proposed an experiment inhibiting Rabi oscillation, 
not relaxation, by frequent measurements. 

We emphasize that our experiments really demonstrate the suppression of relaxation, al- 
beit artificial, and are not mere suspension of Rabi oscillation. 19 ' 20 Although our techniques 
employed in this work are based on the well-known spin-decoupling technique in NMR, 18 we 
propose the new usage of this technique in understanding relaxation phenomena experimen- 
tally and in providing a test bench to develop indispensable techniques in quantum information 
processing. 

We list preceding experimental works dealing with engineered noise here. Kohmoto et al. 22 
analyzed spin relaxation induced by experimentally generated classical random field. Viola et 
al. 23 demonstrated the noiseless subsystem with NMR, in which collective noise was engineered 
through gradient-diffusion method. Kwiat et al. 2A employed the collective artificial noise in 
their decoherence-free subspace experiments with linear optics. Kielpinski et al. 25 applied 
collective noise to ions by irradiating laser light on the ions. Since our noise is generated 
through interaction between qubits, our experiment is closer to realistic situation than theirs. 
4-2 Two-qubit system 

In the rest of this paper we shall study a two-qubit system. Each qubit is referred to 
as qubit 1 and qubit 2, respectively. Qubit 1 is regarded as a system while qubit 2 as an 
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environment coupled to qubit 1. We use the Hamiltonian 



H = JI Z <g) 1, 



Zl 



(35) 



where 1^ = Okj2 (k = x, y, z) and J specifies the strength of the interaction between the two 
qubits. We assume J > without loss of generality. The Hamiltonian of a two-spin molecule 
in a proper rotating frame has this form as will be discussed in Section 5.1. 
When the Hamiltonian (35) acts on states <8> |0) and \tp) <8> |1), it yields 



respectively. Thus the Hamiltonian H describes an effective magnetic field acting on qubit 1. 
The effective magnetic field is J/2 when qubit 2 is in the state |1) and it is — J/2 when qubit 2 
is in |0). Hence, by flipping qubit 2 randomly, we can realize a random fluctuating field for 
qubit 1, see § 3.2. We will extensively use this fact in the following. 

In the next two subsections, we show two examples of artificially generated phase deco- 
herence phenomena: One is phase decoherence of a single qubit in a transmission line, and 
the other is that in a quantum memory. The difference between these two examples is that of 
the characteristic of the random fields. We emphasize that any phase decoherence phenomena 
can be generated by controlling qubit 2. 

4-3 Phase decoherence in a quantum transmission line and its suppression by the bang-bang 



Let us imagine a situation in which a flying qubit passes through a quantum transmission 
line where noise acts on the qubit as it propagates. The noise source is assumed to localize at 
a certain region in the line. 26 

We construct a model which realizes the above situation with the two-qubit system dis- 
cussed in the previous subsection. We regard qubit 1 as a flying qubit and qubit 2 as an 
environment. Suppose that the initial state of qubit 2 is |0). Qubit 2 is flipped to |1) at time 
t±, which corresponds to a position in the transmission line where noise is switched on, and it 
is flipped back to |0) at time t\ + A, which corresponds to a position where noise is switched 
off. We then observe qubit 1 at later time T (> t\ + A). The state of qubit 1 at T is determined 
by applying the phase shift 



on the initial state. Now we regard the time interval A as a random variable. Assume that 
A takes its value in the range < JA < 2ir with uniform probability distribution. Then the 
mixing process (22) yields 



#|Y>)®|0) 
H\i/>)®\1) 



iW) ®io> 

®|i> 



(36) 



control 




(37) 




(38) 
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qubit 1 



m = 2n 



VV^VV^VV^VV^VV^ 



tb 



tb 



tb 



eo 



qubit 2 |0> 



A 
|1> 



time 



|0) 



Fig. 2. Bang-bang pulses, V and V\ are applied repeatedly on qubit 1. Qubit 2 is nipped to |1) at 
ti and is flipped back to |0) at t\ + A. This depicts U\ of Eq. (42) in the case (1). 



Thus, a complete decoherence takes place in qubit 1. 

Let us apply the bang-bang control to qubit 1 in the artificial phase flip channel introduced 
above. Assume that the initial state of qubit 2 is |0). Apply 7r-pulse sequence, V and of 
(23), repeatedly with constant time interval tb ° n qubit 1 during < t < T. Qubit 2 turns 
to |1) at time t\ and turns back to |0) at time t\ + A. The time t\ and the interval A are 
random variables. The value of A varies in the range < J A < 2tt with uniform probability 
distribution. Measure the state of qubit 1 at time T and take an average of the data with 
respect to the random variable A. What is the resulting state ^#(p s ) for this mixing process? 

We consider four cases separately: 

(1) there are even number of 7r-pulses in [0, ti] and even number of 7r-pulses in [t\,t\ + A]; 

(2) even in [0, t\] and odd in [t\,tx + A]; 

(3) odd in [0,ti] and even in [ti,t± + A]; 

(4) odd in [0,*i] and odd in \t\,ti + A]. 

Define eo such that the first 7r-pulse in t\ < t < t\ + A is applied at t = t\ + £q. Similarly, 
define e\ such that the final 7r-pulse in t\ < t < t\ + A is applied at t = t± + A — e\. If the 
number of 7r-pulses in t\ < t < t\ + A is m, 

A = e + (m-l)t h + e 1 . (39) 

Assume that Jt^ <C 2ir so that there are sufficiently many pulses in the interval [tx,ti + A]. 
Under this assumption, the variables eo an d £\ are regarded as random variables taking values 
in the range < e% < tb with uniform probability distribution. 

The time-evolution operator for qubit 1 is calculated for each case as follows. If we put 



LM 



e a 



the time-evolution is generated by (36) as 

e~ iHT \ij) |0) = 



e- iHr \ib) ® II) 



L{-r)\i>) ®|0>, 
L(t)|V) 



(40) 



(41) 
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The sequence of 7r-pulse pairs is represented by alternate insertions of V and of Eq. (23) 
in the time-evolution operator product. For the case (1), the number of 7r-pulses is m = 2n. 
The time-evolution operator is calculated with a help of Fig. 2 as 

U x = V^L(-t h )VL(-(t h - e 1 ))L(e 1 )V^L(t h ) [l(t h )V^ L(t h )v] L(e )L(-(t h - e )) 

= L{t h - {t h -ei) +ei -t h + e - (t h -e )) 

= L(2 ei + 2eo-2t b ), (42) 

where use has been made of the property of V, Eq. (24). For the case (2) with m = 2n + 1, 
we obtain 

U 2 = WL(-(t b -ei))L(ei)r[L(t b )FtL(t b )F] n L(e )L(-(tb-eo)) 

= L(-2 ei + 2eo). (43) 
For the case (3) with m = 2n, 
U 3 = WL(-(tb - ei))L(ei)FL(*b)W [LfojFLfojW] " _1 L(£ )L(-(*b - £o ))FL(-t b ) 
= L(-2ei -2e + 2t h ). (44) 
Finally for the case (4) with m = 2n + 1, 
C/ 4 = yt L (_ tb) y L (_( tb _ £l ))L( £l )W [L(t b )FL(t b )W] n L(e )L(-(t b - e ))VL(-t h ) 

= L(2ei - 2e ). (45) 

By taking average with respect to £o an d £i, and also average over the four cases, the mixing 
process (5) yields 

J Jo ' 

Each term in the parentheses is calculated as; 

L(2e) Ps L\2e) = S(Je)p s S\je) 
I Poo 



(46) 



„iJs „ 

e poi 



while the integral is evaluated as 

1 f th 

tb Jo % Jt\> 
By combining these results we finally obtain 



, vJe = e* Jt » - 1 = sinG7W2) gjJtb/2 _ 



(47) 



(48) 



(49) 



7W2-/- ^ 

Since we have already assumed that Ji b <C 2tt, k should be approximately unity. Comparing 
this result (49) with (38), we see that the bang-bang control suppresses phase decoherence. 



Jt h /2 

Poo «Poi\ 
Jtpw Pu J 



j sin(Jt b /2) ) 2 
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4-4 Phase decoherence in a quantum memory and its suppression by the bang-bang control 

Suppose a qubit sits in a quantum memory device (quantum register). The qubit is exposed 
to a noisy environment and loses its phase coherence. The noise is described by the random 
fluctuating field in the Hamiltonian (18). We assume that the field variable tu(t) in (18) has a 
white noise spectrum. This can also be interpreted as a model of the phase relaxation process 
in NMR. 27 The phase 9 of the qubit evolves in time and is randomly distributed at a later 
time. The distribution function p{9) of 9 is Gaussian, 

^)~^ /2s2 , (51) 



'2txs 

where s 2 is proportional to the evolution time t. 27 

We can also construct a model which realizes the above situation with the previously 
introduced two-qubit system by modifying the random field properly. We regard qubit 1 as 
a qubit in a register and qubit 2 as an environment. Set the initial state of qubit 2 to |0) at 
time to = 0, turn it to |1) at t\, turn it back to |0) at t2, and repeat flipping qubit 2 at t%, t^ 
and so on. Under this manipulation qubit 2 effectively works as a noisy environment acting 
on qubit 1. The time interval between consecutive flippings is denoted as 

A j =t j+1 -t j = A(l + (j = 0,1,2,...). (52) 

Here {£.,•} are independent random variables obeying the probability distribution function 

p {Z J ) = -±=e-tV 2 (53) 



in parallel with (51). A is the average of {Aj}. The parameter a (0 < a < 1/4) characterizes 
variance of the time intervals. To ensure that Aj in (52) is positive, the range of £j should be 
— ^ < £j. However, if a is not too large, the probability of having negative Aj is negligibly 
small. Hence, it is legitimate to extend the range of ^-integration to — oo < £j < oo when we 
take an average. At time i 2 n the evolution operator for qubit 1 becomes 

U 6 = L(-Ai + A 2 + A 2n _! + A 2 „). (54) 

In this case, the mixing process (5) yields 

/oo /'oo poo 

dZl (%2- ■■/ d&n P(6M6)--K6r»)^ft4 ( 55 ) 
-oo J—oo J —oo 

The matrix U^p s ll^ is calculated similarly to Eq. (47). The integral in the off-diagonal com- 
ponents is evaluated as 

/oo poo 
d£ip(£i) / di 2 p(i 2 )e^ J / 2 ^ 1+a ^e^ J l 2 ^ l+a ^ = e -3( JA < (56) 
-oo J—oo 

Thus (55) becomes 

#(ft) /po„ AV.A (57) 

\AVio Pn / 
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at t2n = 2raA, the average of time tin- Then the absolute value of the matrix element poi in 
^#(/9 s ) is multiplied by 

A" = exp{-i(,A a )>! 



exp 



(-ijVA^)=e-WT 2 *. (58) 



The last line defines T 2 * which we call the effective transverse relaxation time since -#(/? s ) 
is calculated only at tin- We note that A is regarded as the correlation time of the artificial 
environment and that the decay of \p§\\ is non-exponential for t ~ A. Thus, we see that phase 
decoherence in the presence of random fluctuating field is characterized by the decay constant 

^ = T^s- m 

It is clear that T 2 * becomes smaller for the larger variance a of fluctuation of the effective 
magnetic field. 

Now let us apply the bang-bang control along with the artificial random fluctuating field. 
Assume that the bang-bang pulse interval % is short enough so that <C a A is satisfied. 
Then the argument of the previous subsection to evaluate ^{p s ) is applicable here as well. 
For t2n = 2nA, there are n random switching of qubit 2 from |0) to |1) and also n random 
switching from |1) back to |0) on average. After one cycle of switching of qubit 2, the element 
poi of the density matrix is multiplied by the factor k given by Eq. (50). Therefore the off- 
diagonal element poi is multiplied by 

K n = K t^/(2A) = e -t^/T* h (6Q) 

at t2n- Equation (60) defines T 2 * b , the decay constant that characterizes the phase decoherence 
under the bang-bang control. It is explicitly given as 



Jtb/2 



sin(Jt b /2) 



-l 



(61) 



Since k approaches 1 from below when % —> 0, T 2 * b approaches oo in this limit. Therefore 
decoherence is suppressed by bang-bang pulses. 

Here we would like to give a remark on the work by Teklemariam et al. 17 They used a 
three-spin molecule as a composite of a one-qubit system and a two-qubit environment. Their 
initial state vector takes the form |x s ) ® \4> e ) ® IV'e)- Their Hamiltonian in our notation is 

H = J 12 I Z ® Iz ® I + Jl3h ® I ® h + J23I ® h ® Iz, 

with which the time-evolution operator U(t) = e~ lHt is defined. They introduced a kick 
operator 

K{Z > ,Q = I®^ Cy ®e Ka \ 
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which acts on the artificial environment. Here £ and £ are random variables, which are inter- 
preted as kick angles. Now the time evolution operator of the whole system is 

■■■K^ 2 ,C2)U(t)K^ 1 ,( 1 )U(t). 

This should be compared with our time-evolution operator (54). Their strategy to manipulate 
the environment is different from ours. A channel associated with their model is defined if Ug ^ 
is substituting into (2). A two-qubit environment is required to simulate an arbitrary one- 
qubit channel. In contrast, we are interested only in the phase decoherence in this paper and 
a one-qubit environment is sufficient for this purpose as was discussed above. It is possible, 
thanks to this simplification, to carry out all the calculations analytically. 

5. Experiments 

We demonstrate generation and suppression of phase decoherence experimentally with an 
NMR quantum computer. A 0.6 ml, 200 mM sample of 13 C-labeled chloroform (Cambridge 
Isotope) in d-6 acetone is employed as a two-qubit molecule. The spin of carbon nucleus in 
chloroform is referred to as spin 1 (qubit 1), while the spin of hydrogen nucleus is referred 
to as spin 2 (qubit 2). Data is taken at room temperature with a JEOL ECA-500 NMR 
spectrometer, 28 whose hydrogen Larmor frequency is approximately 500 MHz. The measured 
spin-spin coupling constant is J/2ir = 215.5 Hz. The transverse relaxation times in a natural 
environment are T 2 ~ 0.30 s for the carbon nucleus and T 2 ~ 7.5 s for the hydrogen nucleus. 
The longitudinal relaxation times are measured to be T\ ~ 20 s for both nuclei. The duration 
of a 7r-pulse for both nuclei is set to 50 /is throughout our experiments. Precision of pulse 
duration control is 100 ns. 
5. 1 Hamiltonian in the rotating frame 

The approximate Hamiltonian of two spins in a heteronucleus molecule, such as Re- 
labeled chloroform (two spins are 13 C and H) in a static magnetic field Bq along the z-axis 
is 

H = -lu 0A I z ® I - u^ 2 I ®I Z + JIz ® h, (62) 

under the secular approximation. 29 Here loq^ = ^Bq, yi being the gyromagnetic ratio of the 
nucleus i, and J is a scalar coupling constant between the spins. The state of the whole system 
evolves following the Schrodinger equation i-^\ip(t)) = H\ip(t)}. If we apply a time-dependent 
unitary transformation 

R = e ~ iuJ °' lht e - iuJ0 > 2ht (63) 
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Therefore, the transformed state \tp(t)) satisfies the Schrodinger equation i^\tjj(t)) = H\ip(t)) 
with the transformed Hamiltonian 




which is nothing but the Hamiltonian (35). The time-dependent operator (63) transforms a 
coordinate system from the laboratory frame to a rotating frame. We will use the rotating 
coordinate system defined with Eq. (63) in the following and the symbol ~ will be omitted 
henceforth to simplify the notations. 

5.2 Demonstration of the phase decoherence in a quantum transmission line 

The model of phase decoherence in a transmission line is implemented with NMR. The 
scheme of the experiment is shown in Fig. 3. Time evolution of spins is depicted in Fig. 3 
(a) in terms of the Bloch vectors viewed from the rotating frame defined with the transfor- 
mation (63). Both spins are set initially in the state |0). This initial state is prepared as a 
so-called pseudopure state. 30 Spin 1 is turned to the x-directon by a 7r/2-pulse along the y-axis 
at t = 0. Spin 2 is flipped by a 7r-pulse at t = t\ and then flipped back to |0) by another 
7r-pulse at t = t\ + A. Spin 1 evolves according to the Hamiltonians (36). Spin 1 precesses with 
the angular velocity —J/2 (J/2) while spin 2 is in the state |0) (|1)). The state of spin 1 at 
t = T is measured via a free induction decay (FID) signal. Figure 3 (b) is a schematic picture 
of the pulse sequence to manipulate these spins. A short bar corresponds to a 7r/2-pulse while 
a long bar to a 7r-pulse. 

The measured quantities via FID signals are the components (a x ,a y ) of the Bloch vector of 
spin 1, which correspond to the off-diagonal elements of the density matrix (13). The measured 



complex amplitudes {a x + ia y = e 1 ^ '-\Ja x + are plotted in Fig. 4. The measurements are 
repeated with variety of spin 2 inversion time A uniformly distributed in the range [0, 2ir/ J] = 
[0, 1/215.5] s. The total number of repetition is 128 in our experiment. 

The open squares □ in Fig. 4 show the measured amplitudes, where only 8 out of 128 am- 
plitudes are shown. Their absolute values are close to unity while their phases are distributed 
in the range [0, 2ir] due to the variation of spin 2 inversion time A. The contribution of the 
trivial phase shift S(— JT/2) in Eq. (37) has been removed when plotting the data. 

To construct a mixing process we take average of measured amplitudes. Each open trian- 
gle A in Fig. 4 is an average of 16 amplitudes. There are 8 = 128/16 open triangles in total. It 
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Fig. 3. Experimental scheme to generate phase decoherence in a flying qubit in a transmission line, 
(a) Motion of the spin Bloch vectors viewed in the rotating frame, (b) The pulse sequence to 
implement the phase flip channel. A short bar denotes a 7r/2-pulse while a long bar denotes a 
7r-pulse. (c) The pulse sequence for the bang-bang control. A regular sequence of 7r-pulses with a 
constant interval ty, is applied on spin 1. 

is found that the absolute values of the averaged amplitudes A become considerably smaller 
than unity. 

The open circle O i n Fig. 4 shows the average of all 128 measured amplitudes. The 
averaged amplitude is close to the origin, which implies vanishing off-diagonal elements of the 
density matrix of spin 1 and therefore is a clear indication of phase decoherence. Thus we see 
that this system works as an artificial phase flip channel for spin 1. 

Next, we apply the bang-bang control to the system qubit. Figure 3 (c) shows the pulse 
sequence to implement the phase flip channel and the bang-bang control. We start applying 
a regular series of short 7r-pulses on spin 1 at t = t\, whose pulse interval is = 0.3 ms. 
The number of 7r-pulses is 16 in each run. The duration of the 7r-pulses is a sum of the pulse 
intervals (0.3 ms) xl5 and the pulse durations (50 /us) xl6 and thus is 5.3 ms, which covers 
the period when spin 2 is in the state |1). The rotation axes of the 7r-pulses for spin 1 are 
cyclically permutated as 

{x,-x,y,-y,-x,x,-y,y) (65) 

in order to reduce influence of possible pulse imperfections. We repeat this procedure 128 
times with randomly varied A's. 

The FID amplitudes in the presence of the bang-bang pulses are indicated by filled 
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Fig. 4. Amplitude of spin 1 measured in the presence of phase decoherence in a transmission line (open 
symbols) and the same with the bang-bang control (filled symbols). □, I : individual amplitude 
measured via FID signals. A, A : amplitude averaged over 16 measurements. 0>* : amplitude 
averaged over 128 measurements (16 x 8 = 128). The rectangle area in (a) is magnified in the 
panel (b). 

squares ■ in Fig. 4. The panel (b) in Fig. 4 is a magnification of a part of the panel (a). 
The absolute values of these amplitudes are almost one. Moreover, their phases are concen- 
trated in a narrow range, \(p\ < 0.25 rad. This is comparable with the theoretical estimate 
Jib = 27r x 215.5 x 0.3 x 10~ 3 ~ 0.4 rad. The filled triangle ▲ in Fig. 4 shows the average 
of 16 amplitudes in the presence of the bang-bang control. The filled circle • in Fig. 4 shows 
the average of whole 128 amplitudes. The average of the whole data is 0.93 + 0.0 i± 0.1 i. The 
error in the y-component originates form the FID phase determination error. 

This magnitude 0.93 may be compared with the theoretical prediction of 0.99. Note that 

we should slightly modify k from Eq. (50) to k = — - — y- - , since the starting times of the 

Jib/2 

bang-bang pulses are fixed at t = t\ in experiments. This discrepancy between our theory 
and experiments should not be taken seriously since the pulse durations (50 /j,s) for qubit 1 
and 2 are finite and are 1/6 of the pulse interval (0.3 ms) of the bang-bang pulses in real 
experiments while they are assumed infinitely short in theory. 

We conclude this section by stating that the bang-bang pulses really suppress decoherence 
generated artificially. 

5.3 Demonstration of the phase decoherence in a quantum memory 

The model of the phase decoherence in a quantum memory is also implemented with 
NMR. We use the same two-spin molecule as that employed in the previous experiment. Both 
spins are initially set to the up-state |0). Spin 1 is turned to the x-axis by a 7r/2-pulse at 
t = to = while spin 2 is flipped from |0) to |1) by a 7r-pulse at t = t^k-i and is flipped back 
from |1) to |0) by a subsequent 7r-pulse at t = t2k {k = 1, 2, 3, • • • ). The rotation axes of these 
7r-pulses are cyclically permutated as in Eq. (65) to reduce undesired influence of imperfections 
in the 7r-pulses. Thus the number of 7r-pulses are set to a multiple of 8. The time intervals 
{Aj = tj + \ — tj} between adjacent 7r-pulses distribute according to the Gaussian distribution 
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(52) and (53). The spin 1 evolves with the Hamiltonians (36). The x- and y-components of 
the Bloch vector of spin 1 at t = T is measured via a FID signal. 

In the first run, depicted in Fig. 5 (a), we put a = and hence the time interval between 
7r-pulses is a constant, Aj = A = 2.0 ms. In this case a regular alternating field acts on spin 1. 
If the 7r-pulses and the spin dynamics were perfect, there would be no decoherence. How- 
ever, in reality, it is impossible to avoid pulse imperfection, measurement errors and intrinsic 
decoherence. Figure 6 shows that decoherence takes place even in the system under regular 
pulses. The measurement with a = is necessary as a reference to the other measurements 
with We may claim that decoherence is enhanced by the random fluctuating field if we 

observe faster decoherence in the measurement with a/0 than that with a = 0. The magni- 
tudes of measured amplitudes under the pulses with a = are plotted as filled squares ■ in 
Fig. 6. The measured decoherence factor for a = is 

e -T/T* = o_45 at T = 100 ms. (66) 

Next, we modulate the time interval of spin 2 flippings randomly. The corresponding 
pulse sequence is shown in Fig. 5 (b). The variance of the intervals in Eq. (52) is chosen 
to be a = 0.10, 0.15, 0.20 and 0.25. The average of the intervals is A = 2.0 ms. A series 
of random variables 3 = (£oj£i>""" is prepared according to the Gaussian distribution 

(53) . Then a series of intervals (Ao, Ai, • • • , A r ) is defined via (52) and then the amplitude 
a x + ia y of spin 1 is measured at time t = T. We prepare 128 series {Si, E2, • • • , Hi2s}, repeat 
measurements 128 times and take an average of 128 measured amplitudes for each values of 



(a) spin 1 
spin 2 

(b) spin 1 
spin 2 

(c) spin 1 
spin 2 



time 



r r — - 1 — - 1 1 — — f 



time 



time 



Fig. 5. Pulse sequences to implement artificial phase decoherence in a quantum memory. Both spins 
are initially in the up-states. A short bar is a 7r/2-pulse while a long bar is a 7r-pulse. (a) A regular 
series of pulses with a fixed interval A = 2.0 ms is applied on spin 2 as a reference to other 
experiments, (b) A series of pulses with random intervals is applied on spin 2. The variance a of 
pulse intervals is adjustable, (c) The bang-bang control pulses are applied on spin 1 while random 
pulses are applied on spin 2. 
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a and T. The magnitude of the averaged amplitude is plotted as a function of T in Fig. 6. 
The correspondence between the symbol and the variance parameter a is (v : a = 0.10), 
(A : a = 0.15), (O : ot = 0.20), (□ : a = 0.25). The plotted data show exponential decrease 
in the magnitude of the amplitude a x +ia y . It is clearly seen that a larger variance a introduces 
faster decoherence in spin 1. The decoherence factors are read from the slopes of the graphs 
in Fig. 6 as 

e -T/T* _ 0.31,0.20,0.10,0.05 at T = 100ms 
for a = 0.10, 0.15, 0.20, 0.25. (67) 



Let us introduce a dimensionless quantity 

R{a) = -\\J e ~ T ^}^ \ (68) 
a 2 J /-'a| a=0 J 

to compare the measured data with the theoretical estimation. The theoretical prediction (58) 
yields 



1 



^theory = ~J 2 AT = 46 (69) 
o 

for J = 2ir x 215.5 Hz, A = 2 ms, T = 100 ms, independently of a. On the other hand, by 
substituting the measured values (66) and (67) into (68), we obtain 

R(a) = 37, 36, 38, 35 
for a = 0.10, 0.15, 0.20, 0.25. 




Fig. 6. Time evolution of the magnitude \a x + ia y \ of spin 1. Decrease in these magnitude indicates 
phase decoherence. The correspondence between the symbol and the variance a of the pulse interval 
distribution is (■ : a = 0), (y : a = 0.10), (A : a = 0.15), (Q : a = 0.20) and (□ : a = 0.25). 
The broken line is the least square fit for {■} while the solid line is that for {□}. 
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We observe that the value R(a) is almost independent of a, implying that the decoherence 
rate (T 2 *) _1 is proportional to a 2 as predicted in (59). 

We conducted a numerical simulation and found that calibration error and spatial inho- 
mogeneity of rf pulse fields may lead to apparently longer T 2 than that with perfect 7r-pulses. 
Therefore, we attribute the small quantitative discrepancy between the theory and the exper- 
iments to rf pulse imperfections. It is observed in Fig. 6 that the data points deviate from 
the straight line in particular for the case with a = 0.25 with T > 80 ms. Our numerical 
simulation also shows that fluctuation in averaged amplitude is large in the region where 
the averaged amplitude is small This fluctuation originates from smallness of the statistical 
ensemble, whose size is 128 in our experiment. 

Next, we apply the bang-bang control to spin 1. The pulse sequence to incorporate the 
bang-bang control is shown in Fig. 5 (c). A regular sequence of 7r-pulses with interval ib = 
0.5 ms is applied on spin 1. The rotation axes of these 7r-pulses are cyclically permutated 
as given in Eq. (65). During this run, a sequence of -zr-pulses whose interval fluctuates with 
variance a = 0.25 is also applied on spin 2. We finally measure the amplitude of spin 1 at 
t = T. We repeat the measurement by preparing 128 series of variables {Hi, • • • , Hi2s} as 
input parameters. The magnitude of averaged amplitude is plotted as a function of T with 
the symbol (□) in Fig. 7. The broken line in Fig. 7 is the reference data with regular pulses 
on spin 2 and no pulses on spin 1 (Fig. 5 (a)). The solid line in Fig. 7 is the experimental 
result with random pulses on spin 2 and without the bang-bang control to spin 1 (Fig. 5 (b)). 
Comparing the data points {□} with the solid line, we observe that decoherence is suppressed 




Fig. 7. Magnitudes \a x + ia v \ of spin 1 with the bang-bang control (t^ = 0.5 ms) applied are shown 
with the symbol (□). The broken line and the solid line are the same as in Fig. 6. Series of 7r-pulses 
with variance a = 0.25 is applied to spin 2. 
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by the bang-bang control. The decoherence factor is read from {□} in Fig. 7 as 

e -T/T 2 * = q 35 a t T = 60 ms for t b = 0.50 ms. 



(70) 



The decoherence factor of the reference broken line in Fig. 7 is 



e -T/T* = Qm atT = 60ms fo r a = 0. 



(71) 



Their ratio is 



e -T/T* 



bang 



0.35 



(72) 



e -T/T 2 * 



= 0.56. 



ref 



0.63 



Theoretical prediction (61) leads 



theory — 0.56, 



(73) 



which is in good agreement with experiments, although there are various factors not taken 
into account in our theory. For example, finite pulse durations, pulse calibration errors, and 
adequacy of the assumption of J% <C 2ir are ignored in our theoretical analysis. 

From these results we conclude that the bang-bang pulses have suppressed decoherence 
caused by the interaction with the dynamical environment. 

6. Conclusions 

We have shown in this paper that a two-qubit system can simulate a composite of a 
system (qubit 1) and its environment (qubit 2) so that qubit 1 exhibits phase decoherence, 
provided that the state of qubit 2 is randomized by external manipulation. We have evaluated 
decoherence rates of qubit 1 and also shown that decoherence is suppressed by applying the 
bang-bang control to qubit 1. 

Performing NMR experiments with a two-spin molecule we measured decoherence in a 
clear manner. In the simulation of phase decoherence in a qubit flying through a quantum 
transmission line, our theoretical calculations were consistent with the measured amplitudes. 
Our theoretical calculations qualitatively explained the measured decoherence rates in the 
simulation of phase decoherence in a quantum memory. It was confirmed that the decoherence 
rate (Tjf ) _1 is proportional to the squared variance a 2 of the interval distribution of the pulses 
applied to the environment (qubit 2). In both cases we demonstrated that the bang-bang 
control successfully suppressed decoherence when the interval of successive time reversal 
operations is much shorter than the correlation time of the artificial environment. 

Study of a qubit system as a composite of a system and its environment will help our un- 
derstanding of the mechanism of decoherence and will help further development of techniques 
to suppress decoherence. 
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